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Abstract 

The non-abehan analog of the Landau- Pomeranchuk-Migdal effect is investi- 
gated in perturbative QCD. Extending our previous studies, the suppression 
of induced soft bremsstrahlung due to multiple scatterings of quarks in the 
spinor representation is considered. The effective formation time of gluon 
radiation due to the color interference is shown to depend on the color repre- 
sentation of the emitting parton, and an improved formula for the radiative 
energy loss is derived that interpolates between the factorization and Bethe- 
Heitler limits. 

24.85. +p, 12.38.Mh, 25.75.+r, 13.87.Ce 



Typeset using REVTeX 



* Permanent address. 



1 



I. INTRODUCTION 



Ultrarelativistic heavy ion collisions at future collider energies of BNL Relativistic Heavy 
Ion Collider (RHIC) and CERN Large Hadron Collider (LHC) are expected to reveal novel 
QCD dynamics not accessible at present fixed target energies. When the transverse mo- 
mentum transfer involved in each nucleon-nucleon collision is small, px ^ Aqcd? effective 
models based on meson-exchange and resonance formation are sufficient to describe multiple 
interaction between hadrons. Those interactions lead to collective behavior in low energy 
heavy ion collisions as first observed in Bevalac experiments []T]] and recently at AGS ener- 
gies 1^ . When pt becomes large enough to resolve individual partons inside a nucleon, the 
dynamics is best described on the parton level via perturbative QCD. Though hard parton 
interactions occur at CERN-SPS energies (-E/ab < 200 AGeV), they play a negligible role in 
the global features of heavy ion collisions. However, at collider energies (-Ecm > 100 AGeV) 
the importance of hard or semihard parton scatterings is clearly seen in high-energy pp and 
pp collisions |^. They are therefore also expected to be dominant in heavy ion collisions at 
RHIC and LHC energies [|,|. 

Since hard parton scatterings occur very early (~ 0.01 fm/ c) and their rates are calculable 
via perturbative QCD (pQCD), we proposed in that high pj- parton jets could serve 
as a unique probe of the quark-gluon plasma formed due to copious mini-jet production 
over longer time scales (> 0.1 fm/c). The systematics of jet quenching provides information 
on the stopping power, dE/dz, of high energy partons in dense matter ^j. The stopping 
power is in turn controlled by the color screening scale, fi, in the medium. Thus, jets 
provide information on that interesting dynamical scale in deconfined matter. At very high 
energy densities, that scale is expected to be large compared to confinement scale 0. At 
temperatures, T 3> Tc, for example, we expect fi ~ gT ^ Aqcd- In that case, most partonic 
interactions have high momentum transfers and perturbative QCD methods may apply to 
the calculation of multiple collision amplitudes. 

In Ref. a systematic study of QCD multiple collision theory was initiated with the aim 
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of deriving the non-abelian analog of the Landau- Pomeranchuk-Migdal (LPM) effect. That 
effect, first derived in the case of QED p!l|JT^ long ago, involves the destructive interference 
between Bremsstrahlung radiation amplitudes. It suppresses radiation relative to the Bethe- 
Heitler formula in kinematic regions where the radiation formation time is long compared 
to the mean free path. In QCD a similar effect is also expected because it mainly follows 
from general relativistic uncertainty principle arguments. However, we showed in that 
specific non-abelian effects influence the detailed interference pattern in the case of QCD. 

The LPM effect in QCD is especially important for estimating the energy loss, dE/dz, of 
an energetic parton traversing a dense QCD medium. We note that there exists considerable 
controversy in the literature on the magnitude and energy dependence of the energy loss in 
QCD |T3Hl7| . In this paper we extend our derivation in to clarify further this effect and 
improve our previous estimates of dE/dz. 

In the next section, we consider in detail the problem of induced gluon radiation from a 
spin 1/2 quark suffering two elastic scatterings. This provides an insight into the applicability 
of the effective potential model used in Ref. []TD| to calculate multiple collision amplitudes. 
We show why radiation from the target partons is negligible although such amplitudes 
are absolutely necessary to ensure gauge invariance. On the other hand, gauge invariance 
constrains all gluon propagators to be regulated by the same screening mass in the potential 
model. In section 3, the QCD radiation interference pattern is shown to be expressible as a 
function of an effective formation time that depends on the color representation of the jet 
parton. In section 4, that formation time is used to calculate the average radiative energy 
loss dE/dz and derive a simple formula that interpolates between the factorization and 
Bethe-Heitler limits as a function of a dimensionless ratio depending on the incident parton 
energy E, the screening scale and the mean free path. Finally, section 5 contains a summary 
and closing remarks. 
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II. THE POTENTIAL MODEL AND GAUGE INVARIANCE 



To analyze multiple parton scatterings and the induced gluon radiation, certain simpli- 
fications have to be made for the interaction. Consider the scattering of a high energy jet 
parton in a color neutral quark-gluon plasma. If the average distance, /S.z = A, between two 
successive scatterings is large compared to the color screening length, A ^ fi^^, the effective 
average random color field produced by the target partons can be modeled by a static Debye 
screened potential: 

viA^) = A'^^Me-^''-- = gTU'^^,. (1) 

q "t" /i 

where /i is the color screening mass, T'* are the generators of SU{'i) corresponding to the 
representation of the target parton at x. The initial and final color indices. A, A', of the 
target parton are averaged and summed over when calculating the ensemble averaged cross 
sections. This model potential was used in Ref. |T0[ to calculate cross sections of multiple 



scattering and induced radiation. However, to obtain a gauge invariant amplitude of gluon 
radiation in QCD, every diagram with the fixed number of final gluon lines has to be taken 
into account, including gluon radiation from the target parton line. Radiation from target 
partons cannot be described by the above static model potential. The relative importance 
of different diagrams to the net energy loss depends on the choice of the gauge. In the 
light-cone gauge one expects that radiation from target legs to be negligible as compared to 
that from the high energy beam parton. 

To estimate the importance of target radiation to the energy loss of a fast parton and to 
see how the gauge invariance constrains the potential model, consider the simplest case of 
induced radiation from quark-quark scatterings. The Born amplitude for {pi, ki) {pf, kf) 
through one gluon exchange is, 

XA -^n^T^ T« u{pf)-f^u{p^)u{kf)^^'uiki) 

where A, A', B, and B' are the initial and final color indices of the beam and target partons, 
respectively. The corresponding elastic cross section is. 



daei _ (i)27ra|s_Mr 

dt ~ S2 t2 



(3) 



where C^^j^^ = Cf/2N = 2/9 is the color factor for a single elastic quark-quark scattering 
and s,u, and t are the Mandelstam variables. 

For induced radiation, there are all together three groups of diagrams as shown in Fig. |l|. 
If we rewrite the three amplitudes as 



(4) 



then gauge invariance implies that 



(5) 



1=1 



where e and k are, respectively, the polarization and momentum of the radiated gluon. In 
the soft radiation limit, we can neglect all the terms proportional to fc^ in A4^^^ which does 
not contribute to A4^^^k^{k'^ = 0). The complete and gauge invariant amplitude for induced 
gluon radiation is. 
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where h is the color index of the radiated gluon, and 

Wa',BB' - )^ JSB'^AA' - >^ \AA' 



(6) 



(7) 



are the color matrices associated with gluon radiation from the internal gluon line (Fig. |^). 
The three terms in Eq. correspond to the radiation from the projectile (Fig. |l|a), target 
parton (Fig. [I|b) and the internal gluon line (Fig. |T]c). As one can see, no term alone is 
separately invariant under a gauge transformation ^ + ck^. Only the total amplitude 
is gauge invariant. 



In the potential model, one simply neglects the radiation from the target lines and replace 
the gluon-exchange amplitude associated with the target parton, gT\ji^,u{kf)'~^^u{ki) / {ki — 
kfY, by an effective potential g'^^A'^y^,. The potential given by Eq. ([^) is regularized by the 
color screening mass squared, /i^. However in Eq. (|^), if any of the gluon propagators is 
regulated, gauge invariance is preserved only if all propagators are regulated with the same 
mass squared, /i^. In particular, the two propagators in the internal radiation diagram must 
both be regulated by the same scale. 

While gauge invariance is manifest only if all diagrams are added, we now show that, in 
the y4+ = gauge, only the projectile diagrams contribute significantly to the radiation of 
soft gluons in the dominant kinematic range for the net energy loss. We consider the case 
that a gluon with light-cone momentum and polarization, 

k = [xP+,fcl/a;P+,kJ, 

6= [0,2e-l-k^/xP+,eJ, (8) 
is radiated off a high energy quark with initial momentum 

p, = [P+,0,Oj. (9) 
During the interaction, the beam quark exchanges a momentum 

g = [?^,9",q±] (10) 

with a target quark which has a typical thermal momentum, 

fci = [M,M,Oj, (11) 

with M ~ T in the plasma rest frame. We focus on the limit defined by x ^ 1, however, 
with the condition, xP^ ^ M ^ qj^. By requiring both of the final quarks to emerge on 
mass shell, one finds that, 

^ " -M' 

(q±-k^)^ fcj .^2) 



The final momenta of tlie beam and target partons are, respectively, 

(q_L — k_|_)^ 

Pf = Pi + q - k [{1 - x)P+, ^^p^ , q± - k^], 

kf = k,-q-[M + ql/M, M, -q^]. (13) 

With the above kinematics, one can obtain the momentum elements of the radiation ampli- 
tudes: 

. pi = 2f^; (14) 

K-pf K ■ Pi Kj_ 

k-ki (xP+)2' k-kf xP+M' ^ ' 

In the large xP+ ^ fc^ limit, we see that the magnitude of the matrix elements involving 
target parton radiation in Eq. ( |T5D are much smaller than those involving projectile radiation 
in Eq. (|14D- As we will see below, the LPM effect limits the radiation to x < AyU^/P"^ and 



fcx < Therefore, as long as A/i ^ 1, the contribution to the energy loss in the main 
kinematic range is dominated by the projectile radiation in this gauge. However, as we have 
demonstrated, the small target contributions to the radiation amplitude are amplified if 
is replaced by k^. As a result of this amplification, those amplitudes cannot be neglected 
when considering gauge invariance even though they can be neglected for calculating the 
energy loss. 

Taking into account only the dominant contributions to the radiation amplitude, we have 
the factorized amplitude. 



^AA'-'-BB' 



7^l ~ 2igei 



k_L q± - kj 



T^^,[T^T%B,, (16) 



kl (q±-k^)2_ 

where M.ei is the elastic amplitude as given in Eq. (0), and TZi is defined as the radiation 
amplitude induced by a single scattering. For later convenience, all the color matrix elements 
are included in the definition of the radiation amplitude TZi. With the above approximations. 
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we recover the differential cross section for induced gluon bremsstrahlung by a single collision 



as originally derived by Gunion and Bertsch 

d(7 dcr^i dn^^^ 



dtdyd'^k± dt dyd^k^ 
where the spectrum for the radiated gluon is, 



(17) 



dyd^k^ 2(27r)3C(;)'^'' fci(qx-k^)2- 

In the square modulus of the radiation amplitude, an average and a sum over initial and 
final color indices and polarization are understood. We see that the spectrum has a uniform 
distribution in central rapidity region (small x) which is a well known feature of QCD soft 



(18) 



radiation [jT9|. This feature is consistent with the hadron distributions predicted by Lund 



string models and the "string effects" in e+e^ three jets events p0| , pl| , all being the results of 
interference effects of pQCD radiation. To demonstrate this a little in detail, let us consider 
only the radiation amplitude from the beam quark in Eq. (0), 



7^ 
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^ ' Pi trnarpb\ 



'BB' 
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The corresponding gluon spectrum is (beside a factor 1/2 (27?)^ 



rpa 

-'-A A'- 



(19) 



1 



C. 



(1) 

el 



\n\ 



e-Pi 



e-pf 



2 k ■ Pi k ■ pf^ 



(20) 



_k-pi k-pf_ 

where Cp = [N"^ — 1)/2N and Ca = N are the second order Casimir for quarks in the 
fundamental and for gluons in the adjoint representation, respectively. Note that the first 
term is identical to gluon radiation induced by an abelian gauge interaction like a photon 
exchange and does not contribute to the gluon spectrum in central rapidity region due to 
the destructive interference of the initial and final state radiations. The second term arises 
from the non-abelian interactions with the target partons and is the main contribution to 
the central region. 

Another nonabelian feature in the induced gluon radiation amplitude, Eq. (|TB|), is the 
singularity at = due to induced radiation along the direction of the exchanged gluon. 
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For k± <^ g_L, we note that the induced radiation from a three gluon vertex can be neglected 
as compared to the leading contribution l/k"^- However, at large k±_ ^ q±, this three gluon 
amplitude is important to change the gluon spectrum to a behavior, leading to a 

finite average transverse momentum. Therefore, q± may serve as a cut-off for k± when one 
neglects the amplitude with the three gluon vertices as we will do when we consider induced 
radiation by multiple scatterings in the next section. If one wishes to include the three gluon 
amplitude, then the singularity at = qj_ has to be regularized. As we have discussed 
above, the regularization scheme has to be the same as for the model potential or the gluon 
propagator in daei/dt as required by gauge invariance. In our case, a color screening mass 
/i will be used. 



III. EFFECTIVE FORMATION TIME 



The radiation amplitude induced by multiple scatterings has been discussed in Ref. |]T0| . 
We discuss here the special case of double scatterings to gain further insight into the problem. 
Consider two static potentials separated by a distance L, which is assumed to be much larger 
than the interaction length, In the abelian case, the radiation amplitude associated 

with double scatterings is (see Appendix), 



7^Q^° = ^g 



(21) 



k ■ Pi k ■ p J \k ■ p k ■ pf ^ 

where p = {p'f,Pz, P±) is the four-momentum of the intermediate parton line which is put on 
mass shell by the pole in one of the parton propagators, Xi = (0,xi), and X2 = (^2,^2) are 
the four-coordinates of the two potentials with ^2 = (-22 — Zi)/vz = Lp^/p^. This formula has 
been recently used to discuss the interference effects on photon and dilepton production in 
a quark gluon plasma p2|j23[] . We notice that the amplitude has two distinct contributions 
from each scattering. Especially, the diagram (Fig. |^) with a gluon radiated from the 
intermediate line between the two scatterings contributes both as the final state radiation 
for the first scattering and the initial state radiation for the second scattering. The relative 
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phase factor 

k ■ {x2 — Xi) = uj{l/vz — cos9)L = L/r^k) 

determines the interference between radiations from the two scatterings and is simply the 
ratio of the path length to the formation time defined as 

""^^^ = a;(lM-cos0) ^ 1^- ^^^^ 

The Bethe-Heitler limit is reached when L ^ T^k). In this limit, the intensity of induced 
radiation is additive in the number of scatterings. However, when L <^ T{k), the final state 
radiation amplitude from the first scattering is mostly cancelled by the initial state radiation 
amplitude from the second scattering. The radiation pattern then looks as if the parton has 
only suffered a single scattering from pi to pf. This destructive interference is often referred 
to as the Landau- Pomeranchuk-Migdal (LPM) effect. The corresponding limit is usually 
called the factorization limit. 

The radiation amplitude in QCD is similar to Eq. ( pID , except that one has to include 
different color factors for each diagram in Fig. ^ In the high energy limit, e ■ pi/k ■ pi ~ 
e ■ p/k ■ p c:^ e ■ Pf /k ■ Pf ^ 2e± ■ k^/A;^. The momentum dependence of each contribution can 
be factorized out and the radiation amplitude for diagrams in Fig. |^ is, 

= z2^^^ {(t-[T-,T^])^^, e^'-^ + ([T-,T^]T-)^^, e^'^-} T-^,T-^, , (23) 

where we have included two color matrices from the potentials, and b again represents 
the color index of the radiated gluon. The radiation amplitude from diagrams with three 
gluon vertices has the same phase and color structures as in Eq. (P5|), but the momentum 
dependence cannot be factorized, since each term depends on the transverse momentum 
transfer which differs from one scattering to another. However, since we are interested in 
the soft radiation limit, kj_ < qj_ ^ /i, the contributions from internal gluon line emissions 
can be neglected as shown in ||10|. As we discussed above, however, those amplitudes serve 
to provide an effective cutoff {q±) ~ fi for k±. 
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The extrapolation of Eq. (p3D to the general case of m number of scatterings is straight- 
forward with the result 

= i2g^Tl\^, . ..TZa'^ E (t"™ ■ ■ ■ [T'^NT^] ■ ■ -T-)^^, e^'^', (24) 

The above amplitude contains m terms each having a common momentum dependence in 
the high energy limit, but with different color and phase factors. The above expression is 
also valid for a gluon beam jet, with the corresponding color matrices replaced by those 
of an adjoint representation. In Eq. (|2^, we also assumed that all the potentials have a 
color structure of a fundmental representation. One can generalize to the case in which each 
individual potential could have any arbitrary color representation. However, our following 
results on the gluon spectrum and interference pattern will remain the same. With this in 
mind, we have the spectrum of soft bremsstrahlung associated with multiple scatterings in 
a color neutral ensemble, similar to Eq. ([T8|) , 



"^mp = Cm{k)— —r- , (25) 



where CgJ""* = {Cp/^N)"^ is the color factor for the elastic scattering cross section without 
radiation. Cm{k), defined as the "radiation formation factor" to characterize the interference 
pattern due to multiple scatterings, can be expressed as 



m 



i-1 



(26) 



where the color coefficients, as computed in Ref. [10|, are 



Cii — C^CaN; 

= -^^(1 - ^Y-^-'C^CaN, for J < ^. (27) 

For a gluon beam jet, one can simply change the dimension to that of an SU{N) adjoint 
representation and replaces C2 = Cp hj the corresponding second order Casimir C2 = Ca- 
We then obtain a general form for the radiation formation factor for a high energy parton 
jet, 
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m i—1 



m-r2i?e^^(l-r2) 



i-j-l ik-{xi-Xj) 



(28) 



i=i j=i 



where 



Ca 



N^/{N^ -1) for quarks with C2 = Cp 



(29) 



r2 = 



2C2 



1/2 



for gluons with C2 = Ca 



Similar to the special case of double scatterings, there are a few interesting limits for the 
above general form of radiation formation factor and the induced gluon spectrum. When 
m = 1, Ci{k) = 1. We recover the gluon spectrum induced by a single scattering in 
Eq. ([ISD in the small k± limit. For multiple scatterings in the Bethe-Heitler limit when 
Lij = \zi — Zj\ ^ T{k) for all i > j, the phase factors average to zero and the intensity of 
the radiation is additive in the number of scatterings, i.e., Cm{k) ~ m. In the factorization 
limit, one has Lij <^ T{k) for all i > j. In this case, the phase factors can be set to unity 
and the summations in Eq. (^) can be carried out to give 




In contrast to the Bethe-Heitler limit, the factorization limit is independent of the number 
of collisions as m 00, and the radiation formation factor approaches l/r2 = 2C2/Ca- 
It is interesting to note that the destructive interference for quarks in the fundamental 
representation is so effective that the radiation spectrum induced by many scatterings is 
even slightly less, l/r2 = 8/9, than by a single scattering. For gluon jets, however, the 
interference is not as complete as for quarks. The induced radiation approaches 2 times 
that from a single scattering. Using these values of Cmik) = l/r2 = 2C2/CA in Eqs. ( pSj) 
and (p!8D, the radiation intensity induced by multiple scatterings is proportional to 2C2 as 
compared to Ca in the single scattering case. The gluon intensity radiated by a gluon jet 
is therefore 9/4 higher than that by a quark due to the interference in multiple scatterings. 
This dependence of LPM effect in QCD on the color representation of the beam parton is a 
unique non-abelian effect. As we will discuss in the following, the effective formation time 
of the radiation in a QCD medium should also take this non-abelian effect into account. 





for quarks 



for gluons 



(30) 
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To see analytically how Cm{k) interpolates between the Bethe-Heitler and factorization 
limits, let us average over the interaction points Xj according to a linear kinetic theory. We 



take an eikonal approximation |]10[ for the multiple scatterings so that the transverse phase 
factors can be neglected in the soft radiation limit, k± <^ q± ^ fi. In a linear kinetic theory, 
the longitudinal separation between successive scatterings, Lj = Zj+i —Zi, has a distribution, 

which is controlled by the mean free path. A, of the scatterings. The averaging of the phase 
factors, 



1 



(32) 



l-iX/T{k) 

enables us to complete the summation in Eq. p8|). Neglecting terms proportional to (1— r2)"^ 
for relatively large m, we have, 

1 + X\k) r2[l + x{k)V 
The non-abelian LPM effect in QCD is therefore controlled by the dimensionless ratio of the 
mean free path to an effective formation time, 

X{k) = A/rQCD(A;), (34) 

where the effective formation time in QCD depends on the color representation of the jet 
parton and is related to the usual abelian formation time in Eq.(p2D as, 

(i\ fi\ 2 cosh y 

TQCBik) = r2T[k) = , . (35) 

This formula for the radiation formation factor is illustrated in Fig. ^ as a function of 
T{k)/X for the case of five collisions (m = 5) and shows how Cm{k) interpolates between 
the Bethe-Heitler limit for small value of T{k) / X and the factorization limit for large value 
of T{k)/X. For radiation with average transverse momentum k± fi the additive Bethe- 
Heitler region is limited to rapidities y < log(r2A/i). In general, the radiation formation 
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factor is proportional to m in the limit of large m. We can therefore regard the radiation as 
being additive to the number of scatterings with the radiation from each scattering simply 
suppressed by the factor + X^) due to nonabelian LPM effect. Since this effective 

formation time is a result of the unique color interference effect in QCD, we should use it 
in the following to estimate the radiative energy loss by a high energy parton traversing a 
color neutral quark-gluon plasma. 



IV. RADIATIVE ENERGY LOSS 

We now apply the effective formation time to derive a simple approximate formula for 
the induced radiative energy loss extending our previous result in As shown in the 

previous section the radiation spectrum is given by Eq. (|18D multiplied by the radiation 
formation factor Eq. (pSf). That factor simply restricts the additive kinematic region to 
tqcdI^)/-^ < 1 and the the incremental change in that factor for each successive collision 
can be approximated by a pocket formula dCm{k) / dm ~ ^^(A — TQcoik)) (see Fig. 3). The 
additive radiative energy loss for each collision beyond the first one is then 

/dn 
d^k^dy^^^^k^ cosh y e{X - tqcdIA:)) e{E - k^coshy), (36) 

where rQCD(^) is given by Eq. (^5]), the second 6'-function is for energy conservation, and 
the regularized gluon density distribution induced by a single scattering is, 

drig _ CAOis ql ^^^^ 



d^k±dy 7r2 A;i[(qx -kx)2 + /i2]• 
As discussed earlier, gauge invariance requires that the singularity at kj_ = qj_ in Eq. (p!8D 
must be regularized by the same color screening mass /i as in the elastic cross section in 
the potential model. Since the transverse momentum transfer q± is the result of elastic 
scatterings, we have to average any function /(q±) of q_L by the elastic cross section, 

(/(q±)) = - /f ^9i^/(qx), (38) 
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where s ~ 6ET is the average cm. energy squared for the scattering of a jet parton with 
energy E off the thermal partons at temperature T. For the dominant small angle scattering, 
the elastic cross sections are, 



dql 



I A 1 



(39) 



where Ci = 9/4, 1, 4/9 respectively for gg, gq and qq scatterings. This average can be 
approximated by replacing q'j_ in the numerator of Eq. (^) with its average value, 

SET 



(gi) = In 



(40) 



In the denominator, we simply replace q^^ by fi^ after the angular integration. The remaining 
integration in Eq. (^6D over the restricted phase space approximately leads to the simple 
analytic formula 



rad 



TT 



2^2 



E 



(41) 




In 



-In 



1 + 4^ 



2^2 



1 + 



(42) 
(43) 



In small k± regime, the phase space is mainly restricted by a small effective formation 
time, tqcd < X, which gives the first term proportional to A. For large k±, the radiation 
becomes additive in a restricted phase space constrained by energy conservation. That region 
contributes to the second term which appears to be proportional to the incident energy E. 
However, in the high energy limit the function I2 (xl/ E and hence the radiated energy loss 
grows only as log^ E. 

The above derivation assumed that the mean free path is much larger than the interaction 
range specified by l//i. As we shall discuss below, this is satisfied in a quark gluon plasma 
at least in the weak coupling limit. Therefore, we can neglect the second term in Ji. For 
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a high energy jet parton, E ^ fi, we can also neglect the second term in J2. The resulting 
radiative energy loss reduces in that case to the simple form, 



dErad _ ^Erad (^2^5 / 2 



dz X 71 

which depends on a dimensionless variable 



(44) 



Since we have used gluon spectrum from a single scattering in Eq. (|18D which is valid for 
all values of k±, the full integration over k± results in the logarithmic energy dependence of 
dErad/ dz. This logarithmic dependence was absent in our derivation in Ref. since there 
only the contribution from the very soft k± < fi region was considered. 

We see that the radiative energy loss dE^ad/dz thus obtained interpolates between the 
factorization and Bethe-Heitler limits as a function of the dimensionless ratio ^. In the 
factorization limit, we fix fiX ^ 1 and let E oo, so that ^ 1. In this case, we can 
neglect the second term in Eq. and have. 



dErad C2(ys,2\t ( '^'^'iE\ , . 

^-^(g.)ln(^^j; (46) 

Thus, the radiative energy loss in the factorization limit has only a logarithmic energy 
dependence (in addition to the energy dependence of {q\))- Due to the non-abelian nature 
of the color interference, the resultant energy loss for a gluon (C2 = C a) is 9/4 times larger 
than that for a quark (C2 = Cp)- In the other extreme limit, we fix E and let yuA 00, 
so that ^ ^ 1. In this case, the mean free path exceeds the effective formation time. The 
radiation from each scattering adds up. We then recover the linear dependence of the energy 
loss dErad/ dz ou the incident energy E (modulo logarithms), 

as in the Bethe-Heitler formula. In both cases, the radiative energy loss is proportional 
to the average of the transverse momentum transfer, (g^), which is controlled by the color 
screening mass as in Eq. (|40|). 
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To see more clearly how the factorization limit is approached, we now estimate ^ for 
a parton propagating inside a high temperature quark-gluon plasma. From Eq. (|39| ) and 
the perturbative QCD expressions for the quark and gluon densities at temperature, T, the 
mean free path for 3 quark flavors is 

= ^mP, + ^,9P9 - —r^ ^ '^C(3)^, (48) 

= ^laPi + ^'aaPa ^ ^9 x 7C(3)— , (49) 

where C(3) ~ 1.2. We emphasize that the above mean free path corresponds approximately 
to the color relaxation mean free path, Ac, and not the momentum relaxation mean free 



path, Ap. As shown in Ref. |Q, Ac ~ a^Ap is generally the shorter of the two in QCD. The 
reason why the color relaxation mean free path controls the radiation pattern is that the 
color current responsible for emitting the gluons is coherent only over a distance scale Ac. 
It take a much longer path length to stop a parton. However, unlike in QED, this longer 
momentum relaxation mean free path is irrelevant for nonabelian radiation. 

Using Eqs. (^) and (49) and the perturbative color electric screening mass, /x^ = AirasT'^, 
we see that ^ appearing in the logarithms has a common energy and temperature dependence 
for both quarks and gluons, 

r^E 63C(3) ^ 9 ^ f K a r^n^ 

With the above expression for S,-, we plot the radiative energy loss in Fig. ^ as a function 
of the beam energy inside a plasma at temperature T = 300 MeV with as = 0.3. The 
solid line is the full expression in Eq. (^4]) while the dashed line is the factorization limit 
corresponding to the first term in Eq. (0). We see that Eq. p6|) approximates Eq. (^4]) 
quite well in this parameter range. The energy dependence of the radiative energy loss is 
due to the double logarithmic function in the formula one of which comes from the energy 
dependence of the average transverse momentum (g^) in Eq. (^Ol). 

The energy loss of a quark in dense matter due to elastic scattering was first estimated 
by Bjorken ||2^ and later was studied in detail |^ in terms of finite temperature QCD. 
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For our purpose, a simple estimate taking into account both the thermal average and color 
screening will suffice. In terms of elastic cross sections and the density distributions for 
quarks and gluons in a plasma, we have, 

where v ~ Q'i/Scc' is the energy transfer of the jet parton to a thermal parton with energy 
uj during an elastic scattering, (g^) is the average transverse momentum transfer given by 
Eq. (^0|). Similar to Eqs. (H) and (^91) , we have. 



Pq\ , /P9\ _ 27ra 



-„(S)+-»(S) = ^^^ (52) 



The elastic energy loss of a fast parton inside a quark gluon plasma at temperature T is 
then given by, 

(54, 

For comparison, we plot this elastic energy loss in Fig. ^. In general, it is much smaller than 
the radiative energy loss and has a weaker energy dependence (single logarithmic). 
Using Eqs. (^OD , ( ^61) and ([50|) , the total energy loss can be expressed as, 

dE dEel , dErad Ca^. 2 1 3ET / 9E Svr^a /«x 

It is interesting to note that both the elastic and radiative energy loss have the same color 
coefficient C2. For high energy partons, the radiative energy loss is dominate over the elastic 
one. For E = 30 GeV, T = 300 MeV, and ~ 0.3, the total energy loss for a propagating 
quark is dE/dz ~ 3.6 GeV/fm. Only about 25% of this amount comes from elastic energy 
loss. 
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V. SUMMARY AND DISCUSSIONS 



We extended our previous derivation by considering the role of gauge invariance and 
target radiation in the case of spin 1/2 quarks to improve our estimate of radiative energy 
loss of a fast parton inside a quark gluon plasma. Our main result [Eq. interpolates 
between the factorization and Bethe-Heitler limits, and has unique nonabelian properties. 
The factorization limit is of course consistent with the general bound [T^ imposed by the 
uncertainty principle, but reveals peculiar energy and temperature dependence of the mean 
square radiation transverse momentum controlling that energy loss. The total energy loss 
is very sensitive to the color screening scale in the plasma. The double logarithmic energy 
dependence of dErad/dz is the result of non-abelian aspects of the LPM effect in QCD. The 
same effect should be responsible for the limited gluon equilibration rate as discussed in 



Ref. [|7|. 

Our derivation improves that in |8|JT0| in a number of ways. First, an effective formation 
time tqcd in QCD radiation was used to account for the color interference due to multiple 
scatterings. The dependence of this effective formation time on the color representation of 
the jet parton gave rise to the different color factors, proportional to C2, for the radiative 
energy loss of a quark and gluon. In contrast both are proportional to Ca in the case of a 
single scattering. Secondly, the gluon spectrum including radiation from both the jet line and 
the internal gluon line and regulated consistently with the requirement of gauge invariance 
was used. However, Eq. ( ^4|) is still an idealization to the physically realizable situation 
in nuclear collisions because a number of strong assumptions were made in its derivation. 
The strongest is the extrapolation of pQCD in a regime (7 ~ 1 and the assumption that the 
interaction range is small compared to the color relaxation mean free path. We have therefore 
left with an explicit dependence on /i in dE/dz since strong nonperturbative variations of 
/i(T) occur in the vicinity of Tc. The basic result that dE/dz is proportional to /i^ is however 
very general and consistent with the uncertainty bounds in |jl6|. Therefore, in a separate 



paper ||28|, we will investigate the phenomenological consequences of interesting temperature 
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dependence of fj,{T) suggested by lattice calculations. 
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APPENDIX: 



The radiation amplitude within multiple scattering theory has been derived in Ref. |T0 
in a general from. In this appendix we use the radiation induced by double scattering as an 
example to demonstrate how the general formula arises from the multiple scattering theory. 
We consider the scattering of a high energy particle off potentials as given by Eq. (p. For 
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simplicity, let us first neglect the color indices as in the case of QED. The amplitude for a 
single scattering is then 

M'i> = 2m6{E, - Ef)47rEME,t), (Al) 

^^^^'^^ = ^2!^^'^^^^^^ 4(q)^..fe)e-^''■^ (A2) 

where E = Ei = Ef and the amplitude f{E,t) is defined such that the differential cross 
section is given by 

^ = n\f{E,t)\', (A3) 

and cTj, (Tf are the initial and final polarizations which should be averaged and summed 
over in the calculation of cross sections. One can check that with the definition of A[q) in 
Eq. (|ip, the above formula leads to da/dt = 47raf/(q^ + /i^)^. 

One can similarly write down the amplitude for double scatterings, 

where the energy conservation at each potential vertex sets the energy of the internal line to 
= E = Ef. Amplitudes involving backscattering are suppressed at high energies because 
of the limited momentum transfer that each potential can impart. When the mean free 
path (or distance between two sequential scatterings) is much larger than the interaction 
range (~ l//i) of the potential in Eq. ([^), the singularity in A{q) can be neglected and the 
integration over (with respect to the z direction of X21 = X2 — xi = + r^) gives us 

M^^^ = 2m6{E,-Ef) J ^u^^{pf){-g^)T^,)U^XP^), 
1(2) = 4(P/ -P)7^ 4(P - p,)e-^(P-P>)■^-*(P/-P)■''^ (A5) 



where p = {E, \J E"^ ~ ^i,^±) is the four-momentum of the internal line. One can derive 
the classical Glauber multiple collision cross section from this amplitude by averaging and 
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summing over initial and final state ensemble of the target . In the limit of high energy 
and small angle scattering, one can neglect the phase factor in the above equation and obtain 
the amplitude [as defined in Eq. ([Al[)], 

f2{E,t) = ^J d'b^[-x{h,E)re^^-\ (A6) 

where 1/2! comes from the different ordering of the target potentials and q_L = P/_l — Pi± 
is the total transverse momentum transfer due to the multiple scatterings. The eikonal 
function Xai,a2(}^, E) is defined as the Fourier transform of the single scattering amplitude 
(besides a factor i/2n), 

X.„.,(b,E) ^-J ^e-^^^-^^^u^M 4(q±)«.,(pi). (A7) 

In the definition of the product of eikonal functions, summation over the polarizations of 
the intermediate lines is implied, 

(A8) 

One can generalize the double scattering amplitude to multiple scatterings and sum them 
together to get the total amplitude, 

°° —ir^ 
fiE,t) = Y: fn{E,t) = — / d%-[-x{h,E)re'^--'^ 

n=l 

= — f d\l - e-x(b,i?)]ei<ix b_ ^^g) 
27r J 

This is recognized as the elastic amplitude in the eikonal formalism. One can generalize this 
to the case of Pomeron exchange so that one can obtain both the total and inelastic cross 
sections for hadron-hadron collisions [^,^ . 

For radiation induced by double scatterings, there are contributions from three different 
diagrams as illustrated in Fig. ^. The total amplitude can be written as, 

AlfJ, = 2m6{E, -Ef-uj)J ^uipf)i-g')[Ta + T, + TMP^). (AlO) 
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e.»=E 



4(^ + k-p,)(^ 



g-i(£+k-pj)-xi-i(p/-^)-X2 



Jf 4(p/ + k-^) 



-fc)2+ie f-p+ie 



4f 



- Vi) 



g-i(£-pi)-xi-i(pj-+k-£)-X2 



g-i(^-pi)-xi-i(p/+k-^)-X2 



(All) 



(A12) 



(A13) 



where we have denoted E — Ef — Ei — cu. Similarly as in double elastic scattering, one can 
integrate over and obtain, 

r„ = 4{pf - p)^ 4(p + k - p.) [-^) e-^(p+k-p0.xi-i(p,-p).x2^ (A14) 

r, = 4{pf + k - p')4 4ip' - P^) (g^) e-(p'-P')--^(P/+''-P')-x., (A15) 

where we have used Dirac equations for the spinors and taken the soft radiation limit (ne- 
glecting terms like / Notice that the four momenta for the intermediate line between 
two scatterings before and after the radiation are different. 



p^{E, ^E'-el,e^), 



p' = {E + uj,J{E + uy 



(A16) 



Define 0d = — (p — p^) -Xi — i(p/ — p) •X2 = p^L+^j, -r^-l-pj -Xi — p/ •X2 as the phase factor 

for double elastic scattering, the corresponding phase factors in the above two amplitudes 
for induced radiation become. 



')D + ik ■ xi, 



(A17) 



where we used p'^ pz -\- UJ /vz, Vz — Pz/E and we defined the time components of the four- 
coordinates as ti = 0, ^2 = L/vz- These two terms can be identified with the initial radiation 
of the first scattering and the final state radiation of the second scattering, respectively. 

In the amplitude of the radiation from the middle line, Fb, we can see that there are 
singularities in both of the two propagators. They should both contribute to the integration 
over Iz- To complete the integration, we may use the identity (with k'^ = 0), 
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With £° = E + LJ, the singularity in the second term put the internal line before the radiation 
vertex on mass-shell with four-momentum p'. This contribution corresponds to the initial 
state radiation of the second scattering. For the first term in the above equation, we can 
make a variable change, £' = i — k. With £q = E, the internal line after the radiation 
vertex is put on shell with four-momentum p. This term can be identified as the final state 
radiation of the first scattering. Using identities. 



= 2£ ■ p (A19) 
(/-^)// = 2e-p'/, (A20) 

one has the amplitude of the radiation from the middle line between two scatterings, 

+ MP, + k - p')^ 4ip' - P.) (-|^) (A21) 

Summing all contributions together, we have the total amplitude for radiation induced by 
double scatterings, 

+ r, + r, = - 4(p; _ p)-A ^(p + k - p,) - ^] e^*-+^^-^ 

2ipz \k-pi k-pj 

which has two distinct contributions induced by each scatterings. In the high energy and 
small angle scattering limit, we can neglect the small momentum transfer at each scattering 
as compared to the beam energy E. We also assume the soft radiation limit in which k± <^ 
q±, with q± being the transverse momentum transfer at each elastic scattering. Substituting 
the above expression back into Eq. ( |A1CI|) , the total amplitude factorizes as 

-MS - -Mi?.7^r , (A23) 
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with the double elastic amplitude M^i given by Eq. ( |A5|) and the radiation amplitude 7^^^^ 



by Eq. (|T]). 

In the case of gluon radiation in QCD, we will have different color matrices for different 
diagrams in Fig. |^. With the color indices as defined in Fig. 0, the corresponding color 
factors for F^, Ff,, and Fc are, respectively, 



Tiai rrii 
A. A' J- 



V 



BB' 



( rj-^hrj-^a2rj-^CL-[ 



BB' 



rpai rpa2 
A. Al J- 



(A24) 



When taking the high energy limit, all the contributions have a common momentum depen- 
dence, e ■ Pi/k ■ Pi c::^ e ■ p/k ■ p c:^ e ■ pf/k ■ pf ~ 2e± ■ ]<.±/k'j_. One can immediately arrives at 
Eq. m- 
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FIGURES 



FIG. 1. Diagrams for induced gluon radiation from a single qq scattering. 

FIG. 2. Diagrams for gluon radiation from the quark line induced by double scatterings. 

FIG. 3. The radiation formation factor Cm{k), m = 5, as a function of T{k)/X = 2cosh.{y)/k±X 
for quarks (solid) and gluons (dashed). 

FIG. 4. The energy dependence of energy loss, dE/dz, of a quark with energy E inside a 
quark-gluon plasma at temperature T = 300 MeV. A weak coupling = 0.3 is used. The solid 
line is the full expression and the dashed line is the factorization limit of the radiative energy loss. 
The dot-dashed line is the elastic energy loss. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408344vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408344vl 




Figure 2: 



This figure "figl-3.png" is available in "png" format from: 
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Figure 3 
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